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Abstract 

We present the complete next-to-leading order (NLO) QCD and NLO electroweak (EW) cor¬ 
rections to the lU-boson pair production associated with a hard jet at the LHC including the spin 
correlated VU-boson leptonic decays. The dependence of the leading order (LO) and NLO corrected 
cross sections on the jet transverse momentum cut is investigated. In dealing with the subsequent 
resonant W-boson decays, we adopt the MadSpin method where the spin correlation effect is in¬ 
cluded. We also provide the LO and NLO corrected distributions of the transverse momenta and 
rapidities of final products. The results show that the NLO EW correction is significant in high 
transverse momentum region due to the EW Sudakov effect. 
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1 Introduction 


Both the ATLAS and CMS collaborations at the CERN Large Hadron Collider (LHC) have discovered 
a neutral boson whose properties are compatible with the Standard Model (SM) Higgs boson [Il[2]. 
The next important goal is to understand the nature of the discovered Higgs boson. Meanwhile, to 
make precision measurement of Higgs properties, the study of signals of Higgs-boson production and 
their backgrounds is quite essential. 

The W + jet (V, V = W, Z) productions at the LHC can be used to determine the gauge boson 
couplings and will allow for a better understanding of the electroweak (EW) symmetry breaking. 
Their production cross sections are currently known up to the QCD next-to-leading order (NLO). 
These production processes are basically backgrounds to the precision measurement of the Higgs-boson 
production which is used to search for physics beyond the SM. Since the precision measurements will 
be possible at the LHC Run 2, and the enormous luminosity is permitted to produce gauge bosons 
with large transverse momentum up to TeV scale, the NLO EW corrections to these processes can 
receive large logarithmic contributions which originate from soft/collinear exchange of EW gauge 
bosons in loop diagrams. Therefore, the EW correction could be strongly enhanced by the EW 
Sudakov logarithms. The EW corrections to most of these processes are yet unknown, although they 
are certainly as important as QCD corrections in some colliding channels. 

The NLO QCD corrections to the W~^W~ -|-jet production at the LHC for on-shell W-bosons were 
calculated in Refs. [3H3- While the NLO EW corrections to this process with leptonic W-boson decays 
have not been investigated so far. R is believed that the calculations for the W + jet productions 
with subsequent vector-boson decays including the NLO QCD and NLO EW corrections at the LHC 
are desired, and one can find that this work is listed in the wishlist (Table 3) of Ref. [Bj. 

In this paper, we present the LO and NLO QCD-I-EW corrected integrated cross sections and 
some kinematic distributions for the pp —)• W^W~ -|-jet -|- A process at the LHC, including the W- 
boson leptonic decays by adopting the naive narrow width approximation (NWA) and the MadSpin 
program. The NLO correction can be divided in 0{a^as) EW and 0{a'^al) QCD corrections. The 
paper is organized as follows: In section [21 we provide a general setup of our calculation. In section [31 
we present the numerical results and discussion for the integrated and differential cross sections at the 
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LO and NLO. Finally, a short summary is given in section [H 

2 Calculation strategy 


Since the calculation strategy of the NLO QCD correction to the VF-boson pair production in associ¬ 
ation with a hard jet at the LHC has already been provided in Ref. [3], in this section we give mainly 
the description of the NLO EW calculation. 

2.1 General description 


We adopt the five-flavor scheme and neglect the masses of u-, d-, c-, s-, 6-quark in both the NLO 
QCD and NLO EW calculations [3]. Since the long lifetime of bottom quark induces that its decay 
vertex is away from the primary production point with a resolvable distance, it is possible to exclude 
experimentally these contributions by supposing that the final 6-quark jet can be tagged in principle. 
In other words, by employing 6-tagging technique we can remove the W~^W~ -j- 6-jet events. 

In the evaluation of the NLO QCD corrections in the five-flavor scheme, the singularity from final 
gluon splitting g ^ bb would be incorrectly removed by 6-tagging, which should be originally canceled 
by that from the 6-quark loop in the gluon self-energy. We deal with this problem in the way similar to 
that in Ref. [3]. We consider the two merged bottom quarks with i? < 0.5 as a “light” jet by adopting 
jet algorithm described in next section. If the two bottom quarks are separated as two individual jets 
(i.e., R > 0.5), we abandon this event due to the 6-tagging. 

Since the NLO QCD correction involves the contributions from the real light-quark emission sub¬ 
processes Ui + dj —>• W~^ + W~ + Uk + di and Ui + dj —)• W~^ + W~ +Uk + di (i = 1,2,3 and j,k,l = 1, 2) 
which are related to the CKM matrix elements, we shall keep the non-diagonal Cabibbo-Kobayashi- 
Maskawa (CKM) matrix in the NLO QCD calculation. In this paper, we take the approximation that 
the CKM matrix has block-diagonal form allowing quark mixing only between the first two generations 
since the mixing to the third generation is negligible, i.e.. 


VcKM = Vij ^ 


3x3 


COS 9c sin 9 c 0 
= I — sin 9c cos 9c 0 
0 0 1 


( 2 . 1 ) 


As explained in Ref. [3], the LO contribution for the pp —)■ W~^W -|-jet -|- X process is not influenced 
by the explicit entries of the CKM matrix due to the CKM unitarity. Analogously, the EW one-loop 
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amplitude for 0 —)• W~^W~ + qiqj + g and the tree-level amplitude for 0 —)• W~^W~ + qiqj + g + ^ 
{q = u,d) are dependent on the CKM matrix via only VikVil,VimV*^ and VikV*i^. These combinations 
of CKM matrix elements equal to 6ij after the summation over the intermediate quark flavors due to 
the CKM unitarity. Therefore, the full NLO EW correction to the pp —)• W~^W~ + jet -|- X process 
is independent of the CKM matrix elements, and we set the CKM matrix to the identity in both the 
LO and NLO EW calculations. 

The pp —)• W^W~ + jet -|- X production at the LO involves the following partonic processes, (see 
EiglUfor example) 


qq W+W-g, ( 2 . 2 ) 

gq W+W-q, (2.3) 

gq W+W-q, (2.4) 

bb W+W-g. (2.5) 

In the partonic processes q stands for the light quarks (i.e., q = u,d,c,s). Because 


of the smallness of the (anti)bottom-quark parton distribution function (PDE) in the proton, the 
contributions of the subprocesses with initial 6 -quark, gb W+W-b, gb —>• W+W-b and bb —)■ 
W+W-g, to the hadronic cross section are suppressed with respect to the subprocesses (I2.2l) - ()2.4p . 
As an example, the contribution to the LO total hardronic cross section from the subprocess bb —)• 
W+W-g at the 14 TeV LHC is less than 0.26%. Therefore, it is reasonable to neglect the NLO EW 
corrections from the subprocesses involving initial 6 -quark in our calculation. Since we perform the 
NLO EW calculation for the W+W- -I-jet production at the LHC in a precision up to 0{a^as), the 
lowest order contributions from the partonic processes in photon-quark collision and their NLO QCD 
corrections, which are of the 0{a^) and 0{a^as), respectively, should be included. These partonic 
processes, called as the subprocesses with initial photon, are written as 


-/qiq) W+W g(q), (q = u,d,c,s), (2.6) 

where the contributions from 76 —)• W+W-b and 76 —)• W+W-b are excluded due the 6 -tagging 
technique. 
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Figure 1: The LO Feynman diagrams for the partonic process uu —>• W~^W g. 


In the NLO QCD and NLO EW calculations we apply the FeynArts- 3.7 [7j package to generate 
the Feynman diagrams and their analytic amplitudes in the’t Hooft-Feynman gauge. The reductions 
of algebraic amplitudes are performed by adopting the FormCalc-7.3 program [8]. The scalar and 
tensor one-loop integrals are evaluated numerically by using our developed in-house library based 
on LoopTools-2.8 [8], in which we add an option to switch to quadruple precision codes when the 
numerical instability is detected. We developed also the interface codes to make use of MadSpin 
program PEo] via Les Houches event file mm- 

2.2 EW virtual correction 


The EW virtual correction to the pp —)• W~^W~ -|-jet -|- X process involves the contributions from the 
EW one-loop diagrams for the processes (I2.2p - (I2.4I) and the QCD one-loop diagrams for the processes 
(j2.6l) . which are all at the 0{a^as). In our calculation, both the ultraviolet (UV) and infrared (IR) 
divergences are regularized by adopting the dimensional regularization scheme md = 4—2e dimensions. 
We employ the definitions of the relevant renormalization constants presented in Ref. m- 

We use the on-shell renormalization scheme to renormalize the relevant masses and wave functions 
m- In the Q:(0)-scheme the fine structure constant a is obtained from the coupling of on-shell e-e-y 
three-point function in the Thomson limit. The bare electric charge is related to the renormalized 
electric charge by eo = (1 -f 5Zg)e. By means of the Ward identity and the on-shell conditions, the 
electric charge renormalization constant dZg in the a(0)-scheme can be written as [13j 


= -Uzaa - \tsiiew6ZzA = 

11 1 op^ 


Et^(O) 


— tan 6w 




M| 


(2.7) 


where 9w is the weak mixing angle, and denotes the transverse part of the unrenormalized 

self-energy at four-momentum squared p^. However, in the soft limit p^ -A- 0, the photon self-energy 
induces large logarithmic terms of \og{rin?^/p?) where mj is the light fermion mass and p is the typical 
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scale of a process. In the a(0)-scheme, each external photon leads to a wave-function renormalization 
constant ^5Zaa-, which exactly cancels the large logarithms appearing in the corresponding electric 
charge renormalization constant 5Zf,. While for a LO process whose number of external photons is 
less than the number of EW couplings, the uncanceled large logarithms have to be absorbed into 
the running fine structure constant by adopting the a{Mz)- or G^-scheme. In this work we take the 
G^-scheme, in which ac^ is defined as 


= 


/ M. 


IjZ. w 

vr 


1 - 


‘w 


Ml 


= a(0)(l -|- Ar) 


( 2 . 8 ) 


where G^ is the Fermi constant, and Ar is the NLO EW correction to the muon decay. The explicit 
expression for Ar can be written as m 


Ar = —5Zaa — cot^ 9w 




M| 






. Et^(O) a(0) 
+2cotgvv + 


M'i dvr sin^ 9w 


2sin^%. 


(2.9) 


Then the electric charge renormalization constant in the G^-scheme can be expressed as 


6Ze^ = 6Zf^'> - ^Ar. (2.10) 

From Eas. (|2.7l) . (I2.9p - (l2.10p we may obtain the analytic expression for 5Ze '' in which the large logarith¬ 
mic terms coming from 5Zaa are canceled obviously. After performing the renormalization procedure, 
the UV divergences are canceled exactly. 



Figure 2: Some representative box and pentagon diagrams for the partonic process uu -A W~^W g. 


Some representative box and pentagon Feynman diagrams for the uu -A- W~^W~g partonic process 
are presented in Figl2j The amplitude for each one-loop Feynman diagram can be expressed as a one- 
loop integral. The tensor and scalar 5-point integrals are directly reduced to 4-point integrals m- 
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The A^-point tensor integrals with < 4 are reduced to scalar integrals with the Passarino-Veltman 
algorithm m - The above reduction procedures have already been implemented in the LoopTools- 
2.8 package. In the calculation of the 3- and 4-point tensor integrals, we meet the numerical unstability 
occasionally in the vicinities of the points with zero Gram determinant. To solve this problem in loop 
integration, we use our developed in-house library which can switch to the quadruple precision codes 
automatically when small Gram determinant occurs [mils]. 

2.3 EW real emission correction 

The real photon emission processes are listed below: 

qq —)• W~^W~ (2T1) 

gq{q) W^W~q{q)-i. ( 2 . 12 ) 

In calculating real photon emission partonic processes, e.g., q{pi)q{p 2 ) —^ {p^)W~{PA)g{pb)l{p&), 

we use the two cutoff phase space slicing (TCPSS) method [19] to isolate the soft and collinear IR 
singularities. Two arbitrary small cutoffs dg and 6c are introduced to separate the phase space into 
soft {Eq < i(5s\/I), hard collinear (Eq > ^6sV§, sie or S 26 < <5cs) and hard noncollinear (Eq > 
si6 and S 26 > <5cs) regions, where s is the partonic center-of-mass energy squared and Sij = {pi + pj)'^ 
is the usual Mandelstam invariant. To get IR-safe results at the LO, we should apply the transverse 
momentum cut of pT,j > p^j on the leading jet for W~^W~ -|- jet events. For the real emission 
processes, we merge the two tracks of jets (or jet and photon) into a single track if they are sufficiently 
collinear. This track combination procedure can lead to collinear IR-safe observables for the process 
(|2.12[) . but will induce the soft-gluon divergence problem for the process (12.lip . This problem can be 
solved by adopting an event selection criterion (see section 13.21) in which W^W -|-jet and W~^W 
events are properly dehned. 

2.4 Subsequent IF-boson leptonic decay 

In the LO evaluation of the W~^W~ + jet production followed by IF-boson leptonic decays, e.g., 
pp —>• W^W~ -|-jet ^ + jet -|- X, we employ the following two methods to generate the 

subsequent decays in order to investigate the effect from the resonant kF-boson spin-entangled decays. 
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1. The naive NWA method. In this method, IT-boson is treated as an on-shell particle without spin 
information, and therefore the IT-boson decay products, lepton and neutrino, can be generated 
isotropically in the VF-boson center-of-mass system. The integrated cross section for pp —)• 


W^W + jet —+ jet -|- X can be easily obtained as 


a [pp —)• W~^W + jet —)• + jet + X) 

= a [pp —)• W^W~ + jet + X) X Br {W~ e~i>e) x Br {W~ —)• , (2-13) 


where Br (W i ly) is the branch ratio for the VT-boson leptonic decay W i v. 


2. The Mad Spin method. Since the spin correlation may affect some lepton distributions sig- 
nihcantly, the MadSpin method taking into account the spin correlation was introduced in 
Refs. [HllTO]. The MadSpin program, which is part of MadGraph5_aMC@NLO [20], can be 
used to generate the heavy resonance decay preserving both spin correlation and finite width ef¬ 
fects to a very good accuracy, and is particularly suited for the decay of resonance in production 
events generated at NLO accuracy nnj. 

In the NLO calculation of the pp —>• W~^W~ -|-jet —)> e ~^-|-jet -|- X process, we shall adopt 
only the MadSpin program to deal with the VL-boson leptonic decays, which can keep both spin 
correlation and finite width effects. 

3 Numerical results and discussion 

3.1 Input parameters 

The relevant SM input parameters are taken as |21j 


Mw = 80.385 GeV, Mz = 91.1876 GeV, Mh = 125.7 GeV 


Mt = 173.21 GeV, = 1.1663787 x 10"® GeV^, as[Mz) = 0.118. 


(3.1) 


In the LO and NLO EW calculations we set the CKM matrix to the identity due to the independence 
on the CKM matrix elements as described in last section [3|. While in the NLO QCD calculations we 
shall keep the non-diagonal CKM matrix and take the values of the CKM matrix elements as 



(3.2) 


The TT-boson decay width is obtained as T\y = 2.045294 GeV by using the MadSpin program, and 
the masses of all leptons and light quarks are set to zero. We take the proton-proton center-of-mass 
colliding energy to be 14 TeV and use the NNPDF2.3 QED PDF (NNPDF23-nlo-as-0118-qed set) [22] 
in both the LO and NLO EW calculations. 

3.2 Event selection criterion 


For the W^W~ + jet production, there exist W^W~ + jet -|- 7 four-particle events at the EW NLO. 
In analyzing these four-particle events, we merge the two tracks of photon and jet into a photon/jet 
track when R^j < 0.5, where R.yj = ■\J {y^y — VjY + V'y and yj are the rapidities of photon and jet, 
respectively, and is the azimuthal-angle difference between the photon and jet on the transverse 
plane. In the case that the photon and jet are merged, we define a parameter as = E^/{E^ -|- Ej) to 
classify the prototype three-body events [23l[2l|. If = 0.7, the prototype event is called as a 

W'^W~ + ^ event and rejected, otherwise, it is treated as a W~^W~ +iet event and kept. However, the 
above event selection scheme leads to the uncanceled collinear IR divergence, because the Kinoshita- 
Lee-Nauenberg (KLN) theorem [251126] is not violated only when we integrate over the full z^ region. 
To tackle this problem, the bare quark-photon fragmentation function is introduced to absorb the 
remaining collinear IR divergence. With this method we should subtract the contribution from the 
range of z^y > 2 :^“* , i.e., 







(3.3) 


where the effective quark-photon fragmentation function Siq-^^{z^) can be expressed as |27j 

(^i^) [2,(1 _ 2,)|-. (3.4) 

in the dimensional regularization scheme and TCPSS phase space splitting scheme at the EW NLO. 
The effective quark-photon fragmentation function ^ 5 —>- 7 ( 27 ) describes the probability of a quark 
fragmenting into a jet containing a photon which carries the fraction z^ of the total energy. As shown 
in Ea. ()3.4p . the effective quark-photon fragmentation function consists of two parts, which correspond 
to the perturbative radiation and non-perturbative production of a photon. At the EW NLO, the 
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bare fragmentation function can be written as m 




aQ\ 


1 / 

2 -k er(l — e) Y /Uj 


“1“ /ij) 


(3.5) 


in the MS renormalization scheme. Then we obtain 

aQl 


^9^7 V ^7 


(Zj) — 


27r 


g->'7 v*7 


(z ^)In 


z^ (1 — z^) 6cS 


+ z^ 


+ Dq—^~^(^Z^, fif) 


(3.6) 


where Dq^^{z^, fj,f) is the finite part of the bare non-perturbative quark-photon fragmentation func¬ 
tion. It is parametrized as 


> 7 (^ 7 )/^/) — 


aQq 


.Pg—>7(^7) 111 


P/ 


(1 - z,)Vg 


+ C 


(3.7) 


and has been measured by the ALEPH collaboration [28], where /tq = 0.14 GeV and C = —13.26, and 
Pq^^{zry) is the quark-to-photon splitting function expressed as 




1 + (1 - z^)'^ 


(3.8) 


In our NLO EW calculation, we have checked the independence of the total cross section on the 
two cutoffs 6s and 6c numerically. Eor the subprocesses involving two jets in the hnal state, we merge 
the two hnal jets ji and j 2 into one jet when < 0.5. In analyzing the pp —)• iy'''IT~ -|- jet —)• 

e"'“/x“t'e7'/x+jet-|-A process, we apply the following additional constraints on the hnal jets and leptons: 


\yq\ < 4.5, \ye\ < 2.5, pT,e > 25 GeV, pT,miss > 25 GeV, % > 0.4, Ru > 0.2 , (3.9) 

where pt,£ is the transverse momentum of the hnal charged lepton. As pointed out in Ref. [3], the 
missing transverse momentum is particularly of interest if the W~^W~ -|-jet production is considered 
as a background process for SUSY searches. 

3.3 Total cross sections 

To conhrm the validity of our codes in calculating NLO corrections, we compared our NLO QGD 
results obtained by employing the EeynArts-3.7/FormCalc-7.3/LoopTools-2.8 programs with 
the corresponding ones by applying the MadGraph5_aMC@NLO package in the four-havor scheme, 
and hnd that both are coincident with each other. We also used the same scheme and parameters as 
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[GeV] 

50 100 200 500 1000 

o-LO [fb] 

p = Mw 
^^ = ^^o 

22976(5) 10035(2) 2833.8(6) 208.56(5) 10.971(2) 
20492(4) 8407(2) 2100.1(5) 115.25(3) 4.329(1) 

[fb] 

p = Mw 

= Mo 

30416(67) 13685(50) 4077(15) 348(1) 23.48(8) 

28912(55) 13006(39) 3747(10) 268(1) 13.25(7) 

<5QCD 

p = Mw 

M = Mo 

32.4 36.4 43.9 66.9 114 

41.1 54.7 78.4 132.5 206.1 

'^NLO [fb] 

p = Mw 

M = Mo 

23344(11) 10059(6) 2733(2) 180.8(1) 8.53(1) 
20972(7) 8588(5) 2095(2) 92.3(8) 3.2(2) 

[%] 

p = Mw 

M = Mo 

1.6 0.24 -3.6 -13.3 -22.2 

2.3 2.2 -0.2 -19.9 -26.1 

QCD+EW rr, 1 

^NLO [fb] 

p = Mw 

M = Mo 

30784(68) 13709(51) 3976(15) 320(1) 21.04(8) 

29392(53) 13187(39) 3742(11) 245(1) 12.12(8) 

JQCD+EW [%] 

p = Mw 

M = Mo 

34.0 36.6 40.3 53.5 91.8 

43.4 56.9 78.2 112.6 180.0 


Table 1: The LO, NLO QCD, NLO EW and NLO QCD+EW corrected integrated cross sections and the 
corresponding relative corrections to the W~^W~ + jet production at the 14 TeV LHC by taking /r = Mw and 
/To for some typical values of p™*. 


in Ref. [3] and get the NLO QCD corrected integrated cross section as 31.921(35) pb, which is in good 
agreement with the corresponding result as 31.970(11) pb in Ref. [3]. 

In Tabled] we list the LO, NLO QCD, NLO EW and NLO QCD+EW corrected cross sections and 
the corresponding relative corrections to the W~^W~ + jet production at the 14 TeV LHC separately. 


The relative NLO QCD, NLO EW and NLO QCD+EW corrections are defined as 


QCD 

^QCD ^ ^NLO _ 2 

ctlo 


EW QCD+EW 

^EW ^ ^NLO _ jQCD+EW ^ *^NLO 

<7LO ’ 0-LO 


(3.10) 


The factorization and renormalization scales are chosen as pj = pr = fJ- = Mw and po for comparison, 
where is a dynamic scale defined as po = + YhiPri with i running over the final photon and 

jets. From this table we find that with the increment of the NLO EW correction becomes more 
significant due to the large EW Sudakov logarithms. For = 1000 GeV, the NLO EW correction 
is about —22.2% and —26.1% by taking p = Mw and po, respectively. 

In order to explore the theoretical uncertainty from the factorization and renormalization scales, 
we define the upper and lower relative scale uncertainties as 


V+ = 


maK{a{pf,pr)} 

/ \ -L1 


mm{a{pf,pr)} 

j., Tj— = - - - — 1, 


(3.11) 
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where and run over the following set: 


( 2 , 2 ), ( 2 , 1 / 2 ) , ( 1 / 2 , 2 ) , ( 1 / 2 , 1 / 2 )| 

I V ho ho / J 


(3.12) 


For = 50 GeV, the LO and NLO QCD+EW corrected integrated cross sections at the 14 TeV 
LHC are obtained as clo = 20492^gg^^ fb and respectively [^ . It shows 

obviously that the relative scale uncertainty can be improved by including the NLO QCD and NLO EW 
corrections. In further analysis we set the factorization and renormalization scales as pf = pr = My/. 

3.4 Differential cross sections 

In this subsection we study the differential cross sections for the pp —^ W^W~ + jet e ~^+ 
jet + X process. 

3.4.1 LO distributions 

We calculate the LO kinematic distributions of the final particles after W-boson decays by applying 
both the naive NWA and MadSpin methods as declared in section 12.41 in order to display the spin 
correlation effect. The spin correlation effect in differential cross sections can be acquired by comparing 
the results obtained by using these two methods, and quantitatively described by the relative deviation 


dehned as 



(3.13) 


We present the LO transverse momentum distributions of p and the corresponding relative deviations 


by adopting two event selection schemes in EigsiSKa) and (b) separately. In FiglSKa) we apply only 


the transverse momentum cut of pT,j > 50 GeV on the leading jet. We see from this hgure that 
there exists discrepancy between the two PT,fi- distributions obtained by adopting the MadSpin 
and the naive NWA methods separately. The discrepancy is particularly obvious in the low Pt,^- 
{PT,fj.- < 30 GeV) and high > 100 GeV) regions, and the relative deviation for the Pt,^i- 

distribution can reach 19.8% in the plotted Pt,^- range. In FiglSKb) we take not only the leading jet 
transverse momentum cut of pxp > 50 GeV, but also the additional constraints in Ea. (|3.9p in the 
event selection, and we see that the relative deviation varies from —7.8% to 10.1% in the plotted 
region. 

^An integrated cross section with upper and lower scale uncertainties is denoted as (rifJ-o, Po)^^- 
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Figure 3: The LO transverse momentum distributions of pT and the corresponding relative deviations by 
applying the naive NWA and MadSpin methods with (a) pT,j > 50 GeV, and (b) prj > 50 GeV and the 
additional constraints in (13.91) . 


In FigsjU^a,b) we depict the LO transverse momentum distributions of e"*" and the corresponding 
relative deviations by applying the naive NWA and MadSpin methods. The results in FigUl^a) are 
obtained with only the constraint ofppj > 50 GeV on the leading jet. It shows that the spin correlation 
effect exists distinctly in the PTe+ distribution, and the relative deviation varies from —17.3% to 24.9% 
in the plotted PT,e+ region. In Figj3{b) both the cut of ppj > 50 GeV and the additional constraints 
in Ea. (|3.9p are taken in the event selection. We can see from the figure that the relative deviation 
between the MadSpin and the naive NWA predictions changes from —14.8% to 16.3% in the PT,e+ 
range of [25, 150] GeV, and reaches its maximum at the position of PT,e+ = 25 GeV. 

Figsl5][a,b) demonstrate the LO distributions of the azimuthal-angle difference (l)e+^- between 
the hnal and on the transverse plane and the corresponding relative deviations by adopting 
the naive NWA and MadSpin methods. Figsi5][a) and (b) are obtained by taking the same event 
selection criteria as in FigsjU^a) and (b), respectively. We find that the distribution from the 

MadSpin method has sizable deviation compared with that from the naive NWA method, and the 
relative deviation increases with the increment of 4>e+p.-- Figl5][a) the relative deviation between 
the MadSpin and the NWA predictions is in the range of [—6.0%, 6.4%]. Fig(5][b) demonstrates that 
varies from —6.6% to 3.1% in the plotted range, and its absolute value can exceed 5% in the 
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PT,e- (GeV) 


Figure 4: The LO transverse momentum distributions of and the corresponding relative deviations by 
applying the naive NWA and MadSpin methods with (a) pT,j > 50 GeV, and (b) prj > 50 GeV and the 
additional constraints in (13.91) . 


region of (pe+p.- < 1-5- 




Figure 5: The LO distributions of the azimuthal-angle difference 4>e+p- between the final e"*" and fj,~ and the 
corresponding relative deviations by applying the naive NWA and MadSpin methods with (a) pxj > 50 GeV, 
and (b) pxj > 50 GeV and the additional constraints in (13.91) . 


In Figsi3Ka,b), FigslD[a,b) and Figsl5]^a,b) we see clearly that the spin correlation effect is demon¬ 
strated obviously in the LO PT , p-i PT , e + 4 > e + p - distributions, and the MadSpin program is an 
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efficient tool in handling the spin-entangled decays of resonant 14^-boson in an efficient and accurate 
way. In further investigations at the QCD NLO and EW NLO, we adopt only the MadSpin method 
in event generation for the pp —)• W^W~ -|-jet —)• + jet -|- X process. 


3.4.2 NLO QCD and NLO EW corrected distributions 

In this subsection, we discuss the NLO QCD and NLO EW corrected differential cross sections for 
the pp —)■ W~^W~ -|- jet -|- X process with the transverse momentum constraint of pT,j > 50 GeV 
on the leading jet. Eor the distributions of the subsequent lE-boson decay products, the additional 
constraints in Ea. ()3.9p are also applied in the event selection. The relative NLO QCD and NLO EW 
corrections to the differential cross section da/dx are defined as 



(3.14) 


Figsl6](a) and (b) describe the LO, NLO QCD and NLO EW corrected transverse momentum 
distributions of and W~ , respectively. The relative QCD corrections to the Pt,w+ Pt,w- 
distributions decrease from 45.2% to 14.8% and from 41.4% to 18.3% with the increment of Ptw^ 
from 0 to 300 GeV, respectively. The relative EW corrections to the Pt,w± distributions decrease 
with the increment of Pt,w^: turn to be negative when Ptw^ > 180 GeV. Due to the large EW 

Sudakov logarithms, the EW correction becomes significant in high Ptw^ region. The relative EW 
corrections to the Pt,w+ Pt,w- distributions are about —4.9% and —3.6%, respectively, at the 
position of Pt,w± — 300 GeV. 

In Figsl3(a) and (b) we present the LO, NLO QCD and NLO EW corrected rapidity distributions 
of W~^ and W~ separately. We can see that the rapidity distributions of IV“-boson exhibit similar 
behavior with W~^. The relative QCD corrections to the yw+ and Pw- distributions increase from 
27.6% to 42.2% and from 28.3% to 41.9%, with the increment of \yw±\ from 0 to 3. For both yy/+ 
and y-!^- distributions, the relative EW corrections are less than 3% in the plotted rapidity range. 

In FiglH]we depict the LO, NLO QCD and NLO EW corrected W-pair invariant mass distributions 
and the corresponding relative corrections. With the increment of Mww from 180 to 800 GeV, the 
relative QCD correction decreases from 39.2% to 7.4%. The relative EW correction is always positive in 
the plotted range and can reach the maximal value of about 5% due to the large positive contributions 
of the subprocesses with initial photon. 
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Figure 6: The LO, NLO QCD and NLO EW corrected transverse momentum distributions of (a) W~^ and (b) 
W-. 
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Figure 7: The LO, NLO QCD and NLO EW corrected rapidity distributions of (a) W'^ and (b) W . 

The LO, NLO QCD and NLO EW corrected transverse momentum distributions of ii~ and 6"*“ are 
shown in Figsl^^a) and (b) separately. From the two figures we see that the relative QCD corrections to 
thepr^^- and PT,e+ distributions vary in the ranges of [22.3%, 33.4%] and [9.3%, 34.8%], respectively, 
as PT,fj.- and PT,e+ range from 25 to 150 GeV. At the position of pT = 150 GeV, the relative EW 
corrections to the PT,fi- and Pr,e+ distributions are about —3.5% and —1.0%, respectively. 
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Figure 8: The LO, NLO QCD and NLO EW corrected VF-pair invariant mass distributions. 
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Figure 9; The LO, NLO QCD and NLO EW corrected transverse momentum distributions of (a) ^ and (b) 

e + 


In Fig JlOl we present the LO, NLO QCD and NLO EW corrected distributions of the missing 
transverse momentum. We find that the relative QCD correction to the pT,miss distribution is stable 
and varies form 25.6% to 32.7% for pT,miss S [25, 150] GeV. The relative EW correction to the 
Pr,miss distribution varies in the range of [—2.3%, 2.7%], and reach its maximum at the position of 
Pt, miss ~ 90 GeV. 
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Figure 10: The LO, NLO QCD and NLO EW corrected distributions of the missing transverse momentum. 

4 Summary 


In this work we present the full NLO QCD and NLO EW corrections to the W~^W~ + jet production 
at the 14 TeV LHC followed by the lE-boson leptonic decays. Our results show that the NLO EW 
correction is significant in high energy region due to the EW Sudakov effect which can be most probably 
detected in the LHC experiments. Eor example, with the jet transverse momentum constraint of pTj > 
1 TeV on the leading jet, the NLO EW correction to the production rate for pp —)• W~^W~ + jet + V 
can reach about —22.2% by taking p = Mw- In the calculations of the pp —)• W^W~ + jet —)• 
+jet + V process, we adopt the MadSpin method to include the spin correlation effect 
and find that the results by using the MadSpin program are more accurate than using the naive 
NWA method. Therefore, we conclude that the spin correlation effect should be taken into account in 
precision calculations. We also present the LO, NLO QCD and NLO EW corrected distributions of 
final VE^-bosons and subsequent decay products, and find that the NLO EW correction is significant 
in high lE^-boson transverse momentum and high lE-pair invariant mass regions. 
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